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ABSTRACT 

We study how the matter dispersed when a supermassive black hole tidally 
disrupts a star joins an accretion flow. Combining a relativistic hydrodynamic 
simulation of the stellar disruption with a relativistic hydrodynamics simulation 
of the tidal debris motion, we track such a system until ~ 80% of the stellar mass 
bound to the black hole has settled into an accretion flow. Shocks near the stellar 
pericenter and also near the apocenter of the most tightly-bound debris dissipate 
orbital energy, but only enough to make the characteristic radius comparable to 
the semi-major axis of the most-bound material, not the tidal radius as previ¬ 
ously thought. The outer shocks are caused by post-Newtonian effects, both on 
the stellar orbit during its disruption and on the tidal forces. Accumulation of 
mass into the accretion flow is non-monotonic and slow, requiring ~ 3 — 10 x the 
orbital period of the most tightly-bound tidal streams, while the inflow time for 
most of the mass may be comparable to or longer than the mass accumulation 
time. Deflection by shocks does, however, remove enough angular momentum 
and energy from some mass for it to move inward even before most of the mass 
is accumulated into the accretion flow. Although the accretion rate rises sharply 
and then decays roughly as a power-law, its maximum is ~ O.lx the previous 
expectation, and the duration of the peak is ~ 5x longer than previously pre¬ 
dicted. The geometric mean of the black hole mass and stellar mass inferred 
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from a measured event timescale is therefore ~ 0.2x the value given by classical 
theory. 


Introduction 


When a star comes close enough to a black hole, strong tidal forces can tear off a 
significant amount of mass from the star, and possibly tear apart the star completely. These 
“Tidal Disruption Events” (TDEs) are thought to be common in the centers of galaxies 
(Komossa & Bade 1999 Wang & Merritt 2004) . Although quite uncertain, the estimated 


rate for main sequence stars (MS stars) to be disrupted by supermassive black holes (SMBHs( 
is ~ 10 -5 yr -1 galaxy -1 


)Wang & Merritt 

2004 

Brockamp et al. 

2011 

van Velzen & Farrar 


2014). TDEs are also expected to happen when a white dwarf (WD) passes too close to 


an intermediate mass black hole (IMBH), possibly in a globular cluster (e.g. Krolik & Piran 
2011 Haas et al.|2012 MacLeod et al.|2014 ) or possibly in a galactic center, but the timescale 
is shorter than for a main sequence star disruption by the square root of their mass density 
ratio. 

The closest distance a star in a parabolic trajectory can approach a black hole without 
being disrupted is called the tidal radius R t ~ i?*(M B #/M*) 1//3 , where i?* is the radius of 


the star, M* is the mass of the star, and 1 


the tidal radius can be written (Phinney 


4_bh is the mass of the black hole. More precisely, 
1989) as R t ~ 50 (fc//) 1 / 6 (M*/M 0 ) 2/3 ^ M B ^R g , 


where k is the star’s apsidal motion constant (a function of its internal density profile), / is 
its binding energy in units of and R g = GMbh/c 2 . The ratio k/ f ranges between 

~ 0.02 for fully radiative stars and ~ 0.3 for fully convective stars (Phinne y|l989 ). The index 
£ parameterizes the radius-mass relation: with i?* oc M 1- ^, £ ~ 0.2 for M* < M 0 , but rises 


to ~ 0.4 for higher-mass stars (Kippenhahn & Weigert 1994). For a complete disruption to 


happen, the actual pericenter of the star’s orbit R p must be < R t , and R t must be outside 
the black hole event horizon. This last condition poses an approximate upper limit on the 
black hole mass: for a MS-SMBH encounter, M B h £) 10 8 (/c//) 1 / 4 (M !|< /Mq) 1-3 ^ 2 Mq and 
for a WD-IMBH encounter, M B h £) 10 5 M Q (M*/M 0 ) -1 . In both cases, the actual order- 
unity coefficient of the maximum mass requires a more careful relativistic calculation and 


generically increases by a factor of several for black holes with larger spin parameters (Kesden 


2012 ). 


Consideration of these events by a number of authors (Hills 1975 Frank & Rees 1976 


Lacy et al. 1982 Luminet 1985) led to a concise description of the fate of stellar debris 
in a full disruption as determined by its energy distribution (Rees 1988). The key—and 
plausible—supposition is that the energy per unit mass of the different fluid elements torn 
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off the star has a roughly flat distribution from ~ to ~ +i7 m *, where E m * is the 

binding energy per unit mass, with respect to the black hole, of the matter on the side of 
the star nearest the black hole. The energy distribution is symmetric about zero when, as 
is generally the case, the star approaches on a parabolic orbit. One immediate consequence 
of this energy distribution is that half the star’s mass escapes the potential well of the black 
hole, while the remaining half travels outward on highly elliptical orbits but then returns 
after a considerable delay (a few months for solar-mass main-sequence stars). This model 


predicts that at its peak, the mass-return rate can be very super-Eddington (Loeb & Ulmer 


1997), and the effective temperature associated with the bulk of the luminosity is generally 
in the EUV (Rees|[1988). Once the mass-return rate passes the peak, this energy distribution 


predicts a subsequent mass-return rate that diminishes oc t 5 / 3 (Phinne y||l98 9). For the last 
twenty-five years, a lightcurve oc t~ 5 ^ 3 has been taken as the hallmark of a TDE. 

In fact, however, there are numerous reasons why this conclusion may have been overly 


hasty. One (Lodato et al. 2009 Guillochon & Ramirez-Ruiz 2013b) is that realistic stellar 


structures lead to an energy distribution that isn’t exactly flat, so that the t -5 / 3 variation 
of the mass-return rate becomes valid only at a late stage or possibly not at all during 
the event. Another (Lodato & Rossi 2011) is that even if the bolometric lightcurve is 


oc t~ 5//3 , observations are nearly always done in the optical or near-UV, i.e., on the Raylcigh- 
Jeans portion of the thermal spectrum, where the flux depends only linearly on the effective 
temperature T, not oc T 4 , as does the bolometric flux. The observed lightcurve’s decline 
should therefore be rather shallower than the decline of the bolometric lightcurve. Still 
another is that photon-trapping during the super-Eddington phase of the accretion may 
create a shallower relation between the mass accretion rate and the bolometric luminosity 


(Krolik & Piran 2012). 


In this paper we wish to raise a more fundamental question. Identification of the mass- 
return rate with the mass accretion rate onto the black hole, where the majority of the power 
is generated, demands that the time between the return of a tidal stream and its accretion 
onto the black hole is short compared to the return time of the stream. This is not a trivial 
requirement. Although all the tidal streams initially have specific angular momentum close 
to the initial specific angular momentum of the star, they have orbital energies much larger 
than that associated with a circular orbit having that angular momentum. That is why their 
orbits are highly elliptical: their semi-major axes are all R p (by ratios > {Mbh /M*) 1//3 ). 
To join a disk at a radius r ~ R p requires that the streams’ energy be reduced by this same 


large factor. Rees (1988) argued that relativistic apsidal precession would cause stream 


intersections near R p resulting in strong shocks. As a result of the dissipation of a large 
part of the streams’ orbital energy, the tidal streams would promptly join an accretion disk 
whose outer radius is ~ 2 R p , the inflow time from there to the black hole is less than the 
























4 


characteristic return time. 


However, if the elliptical orbits of the tidal streams are all mutually aligned, coinciding 
only at their percenters, the only place where stream orbits intersect is at the pericenter, 
but the shocks taking place there dissipate a fraction of the orbital energy that is only 


9 2 , where 9 is the angular thickness of the stream orbits about the midplane (Kochanek 


1994 Guillochon et ah 2014). General relativistic precession of the elliptical orbits’ apsidal 


orientation might help (Cannizzo et al. 1990 Hayasaki et ah 2013), but when the star’s 


pericenter is several tens of R g or more, the precession angle is not great enough to engender 
shocks near R p strong enough to dissipate a large fraction of the orbital energy; in addition, 
as we will show here, although apsidal precession does create shocks, they are near the 
apocenter of the orbit, and those shocks produce enough deflection so as to eliminate the 
contribution of apsidal precession to shocks near the pericenter. For similar reasons, Lense- 
Thirring precession, although very likely present, is unlikely to create sufficiently strong 
shocks near the pericenter (Kochanek 1994 Haas et al. 2012 Stone & Loeb 2012). In 


other words, we have no mechanism to realize a key step in the speedy delivery of returning 
matter to the black hole, the quick dissipation of orbital energy. In fact, the situation is 
even worse than this because merely dissipating orbital energy into heat doesn’t solve the 
problem: it means only that the gas is supported by pressure rather than rotation while its 
characteristic distance from the central mass is, to order of magnitude, unchanged. The heat 
must be radiated in order for mass to move inward, and given the large optical depth of a 
star spread over many tens of black hole gravitational radii, that may not be quick. It is 
the goal of this paper to study in detail the hydrodynamics of the returning tidal streams in 
order to identify the locations and mechanisms of orbital energy dissipation, and therefore 
how those streams may ultimately join an accretion flow onto the central black hole. 

Some hydrodynamical simulation studies of these processes have already been carried 
out ( Rosswog et al.||2009 Haas et ah|2012 Guillochon et al.||2014 ), but all previous work has 
been limited to durations a small fraction of that required to track the return of the majority 
of the star’s bound mass. With our technique, it is possible to integrate long enough (~ 12 x 
the orbital period of the most tightly-bound tidal stream) to follow the return of nearly all the 
star’s bound mass. By so doing, we can determine the structure of the system formed upon 
the streams’ return and examine the degree to which it resembles a classical accretion disk. 
Most crucially, we will determine the level of dissipation due to hydrodynamical processes 
suffered by the half of the star’s mass remaining bound to the black hole. The process by 
which the tidal streams form a disk is often called “circularization” in the literature, in the 
expectation that what results is, in fact, a system in which the fluid follows nearly circular 
orbits around the black hole. Another of the goals of this paper is to define this process in 
quantitative terms so that progress toward complete “circularization” can be measured as a 
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function of time. 

To accomplish onr goals, we must solve two distinct problems: the internal dynamics of 
the star as it falls apart under the tidal stresses of the black hole, and the hydrodynamics of 
the tidal streams as they orbit the black hole and interact with one another. We therefore 
divide our computational method into two corresponding sections. In the first, we follow the 
star as its self-gravity is overcome by general relativistic tidal stresses, and internal pressure 
gradients push its matter outward even as some portions are compressed. This calculation is 
conducted in the frame of motion that the center-of-mass of the star would follow if it were 
unaffected by tides; it follows what happens to the star as the star falls toward the black 
hole, swings around, and is tidally destroyed. In the second, we employ a general relativistic 
hydrodynamics code to compute the hydrodynamics of the resulting tidal streams as they 
orbit around the black hole; this calculation is conducted in the black hole frame. The 
initial conditions for the latter are provided from the results of the former by taking the 
distribution of energy and angular momentum ejected from the star across a range of times 
t' and computing where the associated orbits will reach at a single later time t. The details 
of these methods are presented in Section [2] 

For technical reasons explained below, the actual case we treat is a white dwarf disrupted 
by a relatively small (5OOM 0 ) black hole. Fortunately, many of our results can be scaled 
to the much more common case of main sequence stars disrupted by the larger black holes 
generally found in galactic nuclei. 


Simulations 


2.1. Simulation parameters 


In “typical” tidal disruptions, one might expect the star to be a main sequence star of 
~ 1M 0 while the black hole has a mass ~ 10 6 M 0 or more. Unfortunately, when the ratio 
AIbh/M* is this large, the dynamic range of lengthscales and timescales makes simulation 
rather expensive. The ratio between the semimajor axis of the most tightly-bound tidal 
streams and the tidal radius is ~ (ATbh/AT*) 1 / 3 ; the ratio between the gravitational dynamical 
timescales at the two radii is ~ (Mbh/AT*) 1 / 2 . For this reason, we (and a number of previous 
workers in this held) have chosen to study an encounter between a white dwarf and a smaller 
mass black hole. In particular, we have chosen a white dwarf with mass AT* = O.64AT 0 and 
a black hole of mass ATbh = 500AT 0 . Although arguably rarer than main sequence star 


disruptions, white dwarf disruptions are possible and may have been observed (Krolik & 


Piran 


2011 ). 
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We choose the pericenter of the star’s orbit R p to be identical to the tidal radius R t , 
the distance from the black hole at which a star can be fully disrupted. With this definition, 
Rp = Rt ~ (M B h/M*) 1/3 .R*, where R* is the radius of the star. For a white dwarf of the mass 
we have chosen, i?* = 8.62 x 10 8 cm when its internal structure is assumed to be an n = 3/2 
polytrope. For the black hole mass we have selected, we then find R t = 7.94 x 10 9 cm. 
In gravitational units (which we will predominantly use), the fundamental lengthscale is 
R g = GM-qr/c 2 = 7.38 x 10 7 (M B h/ 500M q ) cm, making R t = 108 R g and 77* = 11.7 R g . If 
G — c — 1, R g — M , an abbreviation we will often use. 

Depending on the context, time is best described in any of several units. In the code, 
the time unit is the relativistic one, R g /c = M when G — c — 1. For discussion of tidal 
stream orbits, the characteristic time is the orbital period in the black hole’s gravity for a 
test-particle with an orbital energy equal to (minus) the contrast in gravitational potential 
energy (with respect to the black hole) from the star’s surface to its center when it passes 
pericenter, Ae n = R*R g /R 2 = 1.01 x 10 -3 . This characteristic timescale for our case is 
to = 6.92 x 10 4 M = 170 s and corresponds to the orbital timescale of the most tightly-bound 
matter, whose semi-major axis is ~ 500M. Time measured in these units will be denoted by 
r = t/to- We set the zero-point of time to be the moment at which the star passes through 
pericenter. 

Finally, we choose the spin of the black hole to be 0 for simplicity. With R p ~ 100M, 
the effects of spin would be minor in any case. 


2.2. Stellar disruption and its consequences 


The first step in our calculation is to compute the trajectory of the star as it passes by the 
black hole on a marginally bound orbit with pericenter R p . We do so including all relativistic 
effects. The dynamics of tidal disruption are computed in a local Fermi normal coordinate 
(FNC) frame moving with this trajectory ( Cheng fe Evans|2013 Cheng fe Bogdanovic||2014 ). 
The tidal field is described by a hexadecapole (including terms up to l = 4) multipole 
expansion of the tidal stress (i.e., curvature tensor) in the FNC system. It is then expanded 
to 4th-order in [Rg/Rit)] 1 ^ 2 , where R(t) is the instantaneous distance of the star from the 
black hole, in order to construct a 2PN (second-order post-Newtonian) approximation to 


the relativistic tidal stress. Within a cubical box 877* on a side, we use a version (Cheng & 


Evans 2013) of the intrinsically conservative grid-based code VH1 (Blondin et al. 2012) to 


compute the Newtonian self-gravity of the star and the Newtonian hydrodynamics induced 
by the combination of stellar self-gravity and the 2PN approximation to the relativistic tidal 
stresses. This simulation extends from a time several dynamical timescales before pericenter 
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passage until several dynamical timescales after pericenter passage. 


Throughout this time, we calculate the rate at which mass leaves the box with varying 
amounts of momentum and internal energy. Figure [T| shows the distribution function for 
all ejected mass in the orbital energy-angular momentum plane. The great majority of the 
matter bound to the black hole has specific angular momentum j in the range 0.86-0.93 x 
the specific angular momentum of the star, j 0 ] the unbound matter distribution is similarly 
centered on ~ 1.1 x jo- A small amount of mass has specific orbital energy as small as 
~ — 1.3Ae7v, but there are significant amounts only for e — 1 > — 0.9A6 at (e is the relativistic 
energy per unit rest-mass; e — 1 is therefore the Newtonian specific orbital energy). If this 
bivariate mass distribution is integrated over j, the mass per orbital energy is, as guessed by 


Rees (1988), nearly flat across the range —0.9Ae7v < e — 1 < Ae 


N- 


The fact that specific energy and angular momentum are so well correlated fits well with 
the classical picture that the width of the energy distribution is set by the width in potential 
energy with respect to the black hole. Matter coming from the side of the star nearer the 
black hole is both deeper in the black hole’s potential and has lower angular momentum. 
The contrast in energy per unit mass is, by definition, ~ A~ R*R g /Rp] the associated 
fractional contrast in angular momentum per unit mass is ~ R*/R p , here ~ 0.1, in excellent 
quantitative agreement with our detailed calculation. 


3. Initialization of the global simulation: mass distribution, orbital shapes, 

and simulation details 

One of the unique features in this study is that the initial condition for the global sim¬ 
ulation is constructed using a highly-accurate tidal disruption simulation. Here we describe 
how the information is conveyed from one simulation to the other. 

Throughout the disruption simulation, mass is ejected through the simulation box edges. 
Every 106 M in time (~ 0.1 x the dynamical time at pericenter), we construct 6.4 x 10 5 
pseudo-particles to represent the distribution in momentum and internal energy of the mass 
ejected from the star, and injected into the domain of the global simulation, at that instant. 
Each particle is placed in a cell in (j — jo)/jo x (e — 1)/Aejv space. The mass, momentum, 
kinetic energy, and internal energy of all fluid elements contributing to each of these cells are 
summed; from these, the pseudo-particle is given an entropy in the sense of p/p 1 . Under the 
assumption that the propagation of each pseudo-particle is entirely ballistic, we project its 
location along the geodesic from its injection point to its predicted location at time r = 0.62 
after the star passes pericenter. 





(j - jo )/jo 



0 

5 

0 

5 

0 

5 
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Fig. 1.— Distribution of mass with respect to specific angular momentum (vertical axis) and 
specific orbital energy (horizontal axis), e — 1 is the energy per unit mass after subtracting the 
rest-mass; e — 1 < 0 is bound, e — 1 > 0 is unbound. The color scale for mass is logarithmic. 
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At the end of the local simulation, we determine how many particles carrying how much 
mass, momentum, orbital energy, and entropy are located in each global simulation cell. 
More than 3 x 10' particles are used in the process so that every simulation cell in the path 
of the tidal streams has enough particles to make its average properties well-determined. 
The local rest-mass density in a cell is the total mass of particles it contains divided by its 
volume; the local momentum density is computed similarly from the sum of the particle 
momenta in the cell. Making use of the new density, we convert the mean entropy back into 
internal energy. This information is then used as the initial data for a global, fully general 
relativistic hydrodynamics simulation in the black hole frame. 

The surface mass density resulting from this initial data is shown in Figure [2] In 
the figure, bound gas orbits in a counter-clockwise direction. The coordinates (A", Y) are 
Cartesian coordinates whose origin is at the black hole position and are in units of R g . The 
core of the original star is at the densest point in the plot, (A ,Y) ~ (—800M, — 1700M). 
The radius corresponding to that point (r ~ 2000M) roughly marks the boundary between 
bound debris (inside) and unbound debris (outside). The ratio of the bound to the unbound 
mass is ~ 0.46 : 0.54. 

In order to check the quality of the ballistic transport assumption and our subsequent 
neglect of tidal stream self-gravity, we computed three measures of the tidal streams’ vertical 
structure at early times in the simulation: their actual vertical scaleheight h\ the hydrostatic 
scalehcight due to the black hole’s gravity, /ibh; and the hydrostatic scaleheight due to 
self-gravity h sg . The latter was computed assuming that the gravity due to the gas’s own 
mass was Y rGE, which should be a reasonable estimate when, as is often the case, the 
horizontal width of the stream is > h. We find that h^u S> h sg > h. In other words, tidal 
stream self-gravity near its initial position is much stronger than the vertical component of 
the black hole’s gravity, but, somewhat fortuitously, our streams begin as thin or thinner 
than self-gravity would produce. This structure is transient—the gas’s pressure expands it 
vertically—but the vertical dynamical time when only black hole gravity is present is the 
same as the orbital dynamical time, so the streams travel a significant portion of an orbit 
before they exceed h sg in thickness. By that time, they encounter the shocks described later 
in this paper and become hot enough that vertical self-gravity is no longer significant. 


3.1. Mass Return Rate 

As we mentioned in the Introduction, a central prediction of the classical theory of 
tidal disruption events is that mass should return in the bound tidal streams at a rate that 
depends on time oc f -5 / 3 . We can check the quality of this estimate by following the orbits of 
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Fig. 2.— Initial surface density (in gm cm -2 ) for the global simulation. The black hole 
and the stellar pericenter are at (X, Y) = (0, 0) and (A", Y ) ~ (0,107), respectively. The red 
curve marks the path of the star’s center of mass as it swung counter-clockwise around the 
black hole; the gas orbits in the same sense. 
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our pseudo-particles until all of them have returned to the pericenter at least once. In reality, 
this mass return rate is substantially modified by hydrodynamics, as our simulation explores, 
but it is useful as a standard of comparison. Figure [3] shows the pseudo-particle prediction 
from our simulation of the stellar disruption; at late times it does indeed follow the t ~ 5//3 
power law. However, those “lates times” don’t begin until r ~ 3; by this time, more than 
60% of the bound mass has already returned to the pericenter. Thus, the time-dependence 
of even the mass-return rate can deviate substantially from the usual expectation. 

It is also worth noting that the return rate at the beginning of the global simulation, 
r = 0.62, is negligibly small. This fact supports our choice of this time as the starting point 
for our hydrodynamic calculation of the later evolution of the system. 


Our result differs in several ways from the predictions made on the basis of analytic 


estimates of the specific energy distribution of the tidal streams (Lodato et al. 2009). They 


found that the return rate peaks at what we call to or t = 1; we find the peak is at r ~ 1.5. 
Perhaps more significantly, we find that the rate of return at that peak lies above the rate 
predicted at that time by an extrapolation of the t -5 / 3 dependence at late times, whereas 
it lies below that extrapolation in the earlier analysis. On the other hand, we find (after a 
suitable rescaling to account for the different parameters) rather better agreement with the 


return rate computed by Guillochon & Ramirez-Ruiz (2013a), both in terms of the peak 


time and the shape of the curve; the only contrasts are that our time of peak return rate is 
slightly later than theirs, and the peak return rate we find is a bit higher above the t ~ 5 / 3 


curve. 


3.2. Mass Distribution in Apsidal Angle 


One of the most important properties of the initial condition affecting subsequent global 
evolution of the gas is its spread in apsidal angle of the gas orbits. Figure [4] shows a single 
orbit from one pericenter passage to the next for several pseudo-particles. Unlike the many 
Newtonian models of post-disruption debris evolution (e.g., Guillochon et al. (2014j)), our 
calculation indicates that the orbital apsidal angles for streams with different semi-major 
axes are in general not aligned. There is a swing of ~ 25° between streams created early in 
the encounter (red in the figure) and streams created later (green in the figure), comprising 
the greater part of the mass. Streams leaving the star early tend to have small semi-major 
axes; those leaving later have a broader distribution of semi-major axis. This misalignment of 
orbits induces collisions between the gas streams where their orbits intersect; earlier-arriving 
gas (smaller semi-major axis) passes the pericenter, heads out towards its apocenter, and 
strikes later-arriving gas. Qualitatively similar rotations in apsidal angle were seen in the 
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Fig. 3.— The black curve is the mass-return rate predicted by the pseudo-particle orbits, 
i.e., without allowance for any hydrodynamics. The red curve is oc t ~ 5 / 3 and is normalized 
so that the black curve asymptotes to it. 
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similarly relativistic calculations of Cheng & Bogdanovic (2014), which treated disruptions 
of main sequence stars by larger black holes. 


Figure [5] illustrates the distribution of mass with apsidal angle. We choose a sign 
convention in which negative 0 indicates a tidal stream pericenter before (in the orbital 
sense) the star’s pericenter; positive 0 after the star’s pericenter. In the initial condition, 
the peak of the distribution, at —3°, is moderately narrow (~ 4° FWHM), but the tail to 
positive 0 extends to such large angles (up to +15°) and has such large amplitude that the 
mass fraction in the peak is only ~ 75%, while a full 25% of the mass is in the tail (0 > +1°). 
With ~ 10° contrasts in apsidal direction common, the angle between many possible orbital 
pairs at their near-apocenter intersection is close to 90°. 


The spread in apsidal angles is, in large part, a relativistic effect. Figure [5] also presents 
the distributions resulting from treating different portions of the tidal disruption in Newto¬ 
nian terms. A purely Newtonian treatment results in a distribution whose shape is similar to 
that given by a relativistic treatment, but shifted ~ 3° toward more negative apsidal angles; 
more importantly, it diminishes the mass fraction in the tail toward larger apsidal angles to 
a bit less than 4%. Mixed Newtonian/relativistic treatments give intermediate results. The 
most significant contributor to the Newtonian/relativistic contrast comes from the relativis¬ 
tic apsidal precession of the star’s orbit as it swings through pericenter; the change in apsidal 
direction per pericenter passage is 3ttR g /R p for a highly eccentric orbit, which is ~ 6° for 
our parameters. Although modified by hydrodynamic effects, each later passage through the 
pericenter region results in additional apsidal precession. 


The net result of these post-Newtonian relativistic effects is to create a strong correlation 
between particularly small semi-major axis, i.e., especially tightly-bound and therefore early¬ 
returning matter, and a positive shift in pericenter direction. 


3.3. Global simulation details 


We use the fully conservative general relativistic magnetohydrodynamics code HARM3D 
for our global simulation ( Gammie et al.||2003 Noble et al.||2009 ); there is, however, no mag¬ 
netic field in the simulation reported here. The code uses piecewise parabolic reconstruction, 
second-order Runge-Kutta time stepping, and the Lax-Friedrichs flux to solve local Riemann 
problems. It is therefore second-order accurate in smooth flows. 


We use modified cylindrical coordinates in Schwarzschild spacetime, where the modified 
spatial coordinate variables (aq, aq, aq) are related to ordinary cylindrical coordinate variables 
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Fig. 4.— Black points are a selection of pseudo-particles randomly chosen from the bound 
set to represent gas density in the global simulation’s initial state. Colored ellipses are single¬ 
orbit trajectories for a small sub-sample of the black points; the starting point for each orbit 
is a point of matching color. Red indicates mass that left the star relatively early in the 
disruption; mass that left the star relatively late is shown in green. 
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A0 [deg] 


Fig. 5.— Distribution of mass per apsidal angle in the initial condition. The black curve 
shows the result of our fully relativistic calculation; the blue curve shows what happens if the 
tidal stress is described in terms of Newtonian gravity rather than relativistic but the star’s 
orbit is relativistic; the green curve shows the distribution if the tidal stress is relativistic 
but the star follows a Newtonian orbit; the red curve illustrates what happens in purely 
Newtonian dynamics. Negative </> indicates a tidal stream pericenter before (in the orbital 
sense) the star’s pericenter; positive (j> after the star’s pericenter. 
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(r, z,<t>) by 

r = exp X\ ( 1 ) 

z — ZQSmh.X 2 ( 2 ) 

(p = x 3 . (3) 

The simulation volume encompasses r in = 30 M < r < r out = 5000M, 0 < (j) < 2n, and 
— 1000 M < z < + 1000 M. Outflow boundary conditions are imposed on all the outer 
boundaries of the volume. We do not treat the region inside r in because we expect that in 
real tidal disruptions, MHD effects not treated here will control the system’s evolution this 
close to the black hole, whereas at larger distances hydrodynamics should be more important. 
The outer boundary is chosen so that ~ 80% of the bound mass stays within this radius 
during the simulation. The role of the modified coordinates is to place the grid-cells, 192 
cells spaced at equal intervals on each of the modified coordinate axes, where they are most 
needed within this volume. The radial grid cells have constant Ar/r, but the sinh function in 
the definition of X 2 strongly concentrates the vertical cells toward the midplane. The degree 
of vertical concentration is changed during the course of the simulation: Zo increases from 
0.59 at the start to 2.94 at r — 2, and finally to 21.08 at r = 5. The midplane cell height 
therefore changes from Az = 0.05M to 0 . 2 M to 1 M. Conservation of mass, momentum, 
and energy is maintained during each grid transition. 

The code units of density and internal energy are chosen so that the maximum value 
in the initial condition for each is unity. The unit of all speeds, including sound speeds, 
is c. In these units, our density and internal energy “floor values” are 10 - 16 (r/M ) -3 / 2 and 
10 — 20 (r/M)~ 5 / 2 , respectively. 

With the zero-point of time when the star passes pericenter, the global simulation is 
terminated at r = 12.2 ~ 35min (its actual duration is ~ 11.6 because it doesn’t begin 
until r = 0.62). By the end of the simulation, 81% of the bound mass has returned to the 
pericenter at least once. 

We assume that the gas is adiabatic throughout the simulation. This is an excellent 
approximation for a WD tidal disruption because the cooling time is extremely long: approx¬ 
imately 80 yr, far longer than the simulation duration. However, the adiabatic index changes 
with temperature as radiation pressure becomes more important relative to gas pressure. We 
account for this change by employing the equation of state in the following way. 

Our reasoning begins with the fact that the total pressure p = p gas +Prad = (7 — 1 )u gas + 
u ra d/ 3 where p gas , p ra d , u gaS) u ra d , and 7 are gas pressure, radiation pressure, gas internal 
energy, radiation internal energy, and gas adiabatic index, respectively. We can then write 

p = (T — 1 )u (4) 
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where 


r = 7 


7-4/3 
1 4 ” U gas / U rad 


(5) 


is an “effective adiabatic index” and u = u gas + u rad is the total internal energy. The ratio of 
gas and radiation internal energy u gas /u rad in equation (J5]) is a function of temperature T, 
i.e. u g as = nkT /(y — 1 ) and u ra d = oT l where n, k, and a are the particle number density, 
the Boltzmann constant, and the radiation constant, respectively. One can then derive T for 
a given pair (n, u ) by solving the following equation for T : 


Pgas = nkT = (7 - 1 )u gas = (7 - 1 )(u - u rad ) = (7 - l)(u - aT 4 ) . ( 6 ) 


If 7 , the gas adiabatic index, is set to 5/3, T = 5/3 for u gas » u rad and T = 4/3 for 
u g as « u rad . We find in practice that T is so slowly-varying that it is unnecessary to 
account for its cell-to-cell contrast in the Riemann problem solutions. 


4. Results 
4.1. Overview 

Perhaps the principal result of our simulation is that tidal streams do not quickly and 
easily join an accretion disk at r ~ R p immediately upon their return. Instead, hydrodynamic 
effects only gradually “circularize” their orbits, and the majority of the mass settles at radii 
considerably larger than R p . 

Figures [ 6 ] and [7] show the time evolution of surface mass density and surface internal 
energy density. At early times, the surface mass density largely reflects the paths of the 
tidal streams, but clearly displays a point of intersection near the streams’ apocenter. As 
time goes on, the tidal stream features spread and lose contrast. As they do, the radial 
profile of the surface mass density reaches a slowly-evolving state in which the azimuthally- 
averaged surface mass density is approximately cx r for r < 150M, is roughly flat out to 
r ~ 600M, and then declines steeply at larger radii. Through much of this period, there 
are two prominent ridges of higher density (most noticeable in the panels corresponding to 
5 < r < 11). The origin of these features becomes clearer when the surface mass density 
and surface internal energy density are viewed together. 

The surface density of internal energy is particularly useful for identifying shock lo¬ 
cations, although not all regions of elevated internal energy density contain shocks. For 
example, early in the evolution (r ~ 0 . 8 ), a shock appears close to the stellar pericenter 
(X, Y) ~ (0,100M); this shock can be seen clearly as the bright spot just above the black 
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hole in the r = 2.06 panel of Fig. [7j The result of vertical gravitational focusing (Kochanek 
1994), this shock is often called the “nozzle shock” (Evans & Kochanek 1989 Guillochon 
et al. 2014), we call it “shock 1”. As we will show later, hydrodynamic effects alter the 


structure and even location of shock 1 over time as some material makes multiple passes 


through this region (see Fig. 11 for a pedagogical sketch of this and other shocks’ evolution). 


Shortly after the formation of shock 1, an outer system of shocks forms at large radius 
(r < 1000M) where gas flowing outward after passing through the pericenter region collides 
with newly returning gas. Although visible in the r = 3.51 panel, these shocks are more 
easily seen beginning at r = 4.95: they are the two legs of the feature resembling a “A”. As 
can be seen in Figure [4j in the limit of very narrow stream width, the only shock created in 
a set of perfectly-aligned orbits would be the nozzle shock; outer shocks appear only as the 
result of finite stream width, or, more powerfully, as the result of rotation in apsidal angle. 
Much as in the case of supernova remnants, the encounter between gas on its way out and 
newly returning gas creates a “forward” and a “reverse” shock, which we call, respectively, 


shock 3 and shock 2 (a shock in roughly the same location as shock 2 was noted by Rosswog 


et al. (2009)). Material falling in from large radius encounters the forward shock, shock 3; 
shock 2 acts upon material that has already taken at least one turn around the black hole. 
Initially, shock 2 has only one segment, located very close to shock 3 (see the r = 2.06 and 
r = 3.51 panels of Fig. [7J the faint red feature at large negative Y and small positive A" 
marks these two shocks). However, as gas accumulates behind it, shock 2 splits (see the 
panels for r = 4.95-10.73). One branch stays at roughly the same location and remains 
very close to shock 3, but a new branch moves upstream, swinging ~ 1 radian clockwise in 
azimuthal angle from the stationary branch (see any of the later panels in Fig. [7]). The early 
appearance and, to some degree, the very existence, of shocks 2 and 3 is a consequence of 
the rotation in apsidal angle of the more tightly-bound tidal streams. 

All of the shocks broaden the fluid’s specific angular momentum distribution by deflect¬ 
ing streams away from their ballistic orbits; in these deflections, some streams lose angular 
momentum while others gain it. Those losing angular momentum move inward on more 
eccentric orbits, while those gaining angular momentum traverse less eccentric orbits with 
larger pericenters. Although, as we will discuss in greater detail in Sec. 4.5[ this process 
creates more circular orbits for most of the gas mass, it does not involve significant orbital 
energy dissipation or lead directly to matter joining an accretion disk at radii near R p . In 
fact, after the bulk of the mass has returned (r ~ 4), the surface density peak gradually 
shifts outward from the vicinity of the stellar pericenter to ~ 3 R p ] this shift is due to the 
passage through the inner boundary of matter that has lost angular momentum, so that the 
mean angular momentum in the matter left behind increases. 
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Fig. 6.— Time evolution of surface density from r = 2.06 to the end of the simulation, 
r = 12.12. The time interval between snapshots is ~ 1.5 in these units. The spatial axes are 
in units of M and are centered on the black hole. The surface density scale is logarithmic 
in code units, with a factor of 2000 from greatest surface density to least. One code unit of 
surface density is 1.11 x 10 13 gm cm -2 . Note that the highest surface density locations in 
the r = 2.06 panel can be as much as two orders of magnitude greater than the maximum 
value of the color scale. Such large surface densities are not present at later times. 
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Fig. 7.— Time evolution of internal energy surface density at the same times as the corre¬ 
sponding panels in Fig. [6] The spatial axes are in units of M and are centered on the black 
hole. The internal energy surface density scale is logarithmic in code units, with a factor 
of 1.5 x 10 4 from greatest internal energy surface density to least. One code unit of energy 
surface density is 9.96 x 10 33 erg cm -2 . 
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Removal of angular momentum leading to inflow well within R p can be seen in the 
data of Figure [8j which shows how the radial mass flow F m (r) = — J pv r rdzd(j) depends on 
distance from the black hole for various times. Initially, there is hardly any mass inside 
r ~ R p ~ 100 M and there is therefore no mass flow at small radii; meanwhile, at large 
radii, returning streams rapidly carry mass inward. By r — 3, there is a steady inward flow 
across the range of radii r min = 30 M < r < 125M. This flow is due to the inward deflection 
of streams at shock 1, which initially forms at r ~ R p ~ 100 M. The zone over which the 
inward flow extends gradually stretches outward, as deflections at shocks 2 and 3 also begin 
to contribute; for r > 6, the inflow is in an approximate steady-state out to r < 300M. 

We deliberately do not simulate the dynamics of mass flow at smaller radii than 30M 
because we expect that in this region inflow is controlled by internal MHD stresses rather 
than the shocks and other hydrodynamic processes treated in our simulation. For this reason, 
we distinguish in this paper between mass “accretion” and mass “accumulation”. We reserve 
the former term for mass acquired by the black hole, having suffered enough dissipation of its 
orbital energy to transform ~ 10% of its rest-mass energy into heat and possibly radiation. 
Our simulation therefore treats mass “accumulation”, the capture of mass from the disrupted 
star so that it resides within ~ 100-1000M of the black hole, rather than mass accretion. 

The shocks also dissipate kinetic energy into heat: the increasing importance of pressure 
support causes the geometric arrangement of the gas to grow rounder (see Fig. [6]) and 
thicker from the time the first tidal streams return to the black hole until r ~ 3-6. The 
vertical thickening of the structure can be seen directly in Figure [9j showing the radial 
profile of the azimuthally-averaged disk aspect ratio H/r (H is the scale height of the disk 
= f p\z\dzj f pdz) at several times. The scale height grows rapidly from the beginning of 
the simulation until r ~ 3, and then slowly diminishes for r > 6. The late-time thinning 
is due to the diminishing strength of the shocks combined with adiabatic cooling as the gas 
spreads. Despite the shocks’ gradual weakening, the flow nonetheless remains fairly thick by 
conventional accretion disk standards, with H/r ~ 0.5 for r > 100M and larger at smaller 
radii. 

In the following subsections we treat several of the flow’s chief elements in greater detail. 


4.2. Shock 1 (Nozzle Shock) 

The tidal stream orbital planes have angles with respect to the midplane ~ R*/R p ~ 
(Mbh/Mi.) - 1//3 - If the streams were composed of collisionless particles, they would pass 
through the midplane twice, in most cases once shortly before reaching pericenter and once 
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Fig. 8.— Radial profile of radial mass flow; F m (r ) = — f pv r rd(j)dz. Black, red, green, 
blue, and cyan colors correspond to r =1, 3, 6, 9, and 12, respectively. Positive mass flow 
is inward. The switch in sign between r = 6 and r = 9 at r ~ 800Mbh corresponds to the 
orbital motion of the lump discussed in 


4.3 and also visible in Fig. P9 
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Fig. 9.— Radial profile of azimuthally averaged disk aspect ratio H/r where H is scale 
height defined as f p\z\dzj f pdz. At each radius, H/r is weighted by surface density and 
averaged in azimuthal direction. Black, red, green, blue, and cyan colors correspond to r =1, 
3, 6, 9, and 12, respectively. 
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more shortly after. Because there are streams whose apocenters are both above and below 
the midplane and these streams are actually fluid, this strong convergence has long been ex¬ 
pected to produce shocks near the pericenter (Kochanek||l994). By including high-resolution 
hydrodynamics, we are able to see how pressure forces alter the structures of these shocks 
and prevent the streams from crossing through the midplane. At first, the vertical conver¬ 
gence of the streams merely increases the pressure of gas near the midplane at r ~ R p . Soon, 
however, this high-pressure gas begins to expand vertically; when it does so, it intercepts 
new streams and a shock forms wherever newly-arriving material strikes this high-pressure 
barrier (this process is analogous to the “pancake mechanism” inside a tidally-disrupted star 


described by Brassart & Luminet (2008)). 


The structure of shock 1 can be most clearly discerned in a plot of gas entropy (Fig. 10). 


Region A in this figure is the zone of cold gas arriving at the pericenter region for the first 
time since it left the star. This gas has previously been shocked only in shock 3. As it heads 
toward the midplane from both above and below, it is compressed. After passing through a 
“nozzle” very close to the midplane, it expands on the far side (Region B). Other Region A 
gas, traveling farther from the midplane, encounters the wedge of high pressure at the nozzle 
and is shocked, thus entering Region D. Gas approaching shock 1 at still greater height from 
the midplane (Region C) has higher entropy because this is not its first time through the 
pericenter region, and it has already passed at least once through both shock 1 and shock 2. 
It, too, passes through the shock front to enter Region D. 


The location, size, and strength of shock 1 change over time. In the beginning (r < 1.5), 
shock 1 stays close to the pericenter region for the more tightly-bound streams. Later on, 
shock 1 stretches radially, both inward and outward as the relative contribution of the 
“Region C” gas, the gas that has already been disturbed in either shock 1 or shock 2, grows. 
The trajectories of these streams do not go directly through the pericenter region, as the 
angular momentum distribution of the gas has been broadened by previous shock passages. 
Its vertical extent also changes over time. When it first forms, it extends only a few M 
from the midplane; by r ~ 5, it reaches ~ 50 M above and below the midplane. As shock 1 
stretches, it also becomes weaker: the larger the volume occupied by high pressure gas in 
the nozzle region, the less incoming streams are accelerated toward the midplane by gravity 
and the smoother the overall flow. The shock disappears at r ~ 7, but immediately before 
it does, it spans all the way from the inner boundary of the simulation box r = 30 M to 
r ~ 500 M. 


We close this section with a technical note. Early in the simulation, the gas is still 
extremely cold, making the vertical structure near shock 1 extremely narrow. At r ~ 0.8-1, 
even with our most concentrated vertical cell distribution we have at some times as few as 
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Fig. 10.— Color contours of entropy in the X-Z plane as it cuts through shock 1 at Y = R p 
and r = 2.35. Black arrows show the velocity field of the gas. Different regions are labeled 
by letters and described in detail in the main text. Shock fronts are indicated by black lines. 
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~ 5 cells per scale height in that region. However, the disk rapidly becomes thicker due 
to shock heating (initially mostly at shock 1), and hence we do not think this short time 
period of poor resolution affects the disk evolution. After the disk is thickened enough, we 
deconcentrate the grid at the equator as described in 
nozzle region has H/Az ~ 40 — 50 through nearly all the simulation. 


3.3 With these adjustments, the 


4.3. The outer shocks: shocks 2 and 3 


As described in the previous subsection, shocks 2 and 3 emerge when the most bound, 
and therefore earliest-returning, gas starts to collide with later-arriving gas near the apocen- 
ter of the more tightly-bound gas’s orbit. That the velocities of the two streams are nearly 
perpendicular at the collision point is a consequence of the correlation between apsidal angle 
and tidal stream binding energy. In other words, a combination of two small relativistic ef¬ 
fects (apsidal precession of the stellar trajectory, the exact form of the tidal stress) creates a 
broader spread of orbit orientations, with the contrast largely between the earliest-returning 
streams and those returning later. That early-returning streams are oriented differently from 
the streams arriving soon after promotes the quick emergence of shocks 2 and 3. Stream 
intersections are further promoted by additional apsidal precession when the tidal streams 
return to the pericenter, but fluid effects (notably shock 1) mask that precession. Thus, the 
dissipation in this shock system is a result of small, but consequential, relativistic effects in 
the dynamics of tidal disruption. 


To understand the structure, evolution, and consequences of these shocks, we have 

In this diagram, we designate the newly- 


constructed a schematic shown in Figure 11 


arriving stream gas-1 and the stream that has already swung through the pericenter region 
gas-2. The first step in the shock formation process is a compression of the material in gas-2 
as it decelerates while climbing out of the black hole’s potential well toward apocenter; this 
is a purely kinematic effect and would take place in a collisionless gas in the same way as in a 
fluid. The deceleration creates a density enhancement just prior to the point of intersection 
between gas-1 and gas-2. 


Shocks 2 and 3 form when gas-2 strikes gas-1 from the side. At shock 3, gas-1 is shocked, 
raising its pressure. This high pressure zone creates shock 2, which shocks gas-2. Shocked 
gas-1 and gas-2 then have the same pressure, but are separated by a contact discontinuity 
(left panel of Fig. |Ttj ). Because the density of gas-2 was already elevated by its deceleration, 
its compression in the shock creates a region of even higher density. As time goes on, the 
mass contained in this feature grows, widening the high pressure region between shock 2 
and the contact discontinuity, and moving the location of the shock front farther upstream. 
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Fig. 11.— Schematics of shocks 2 and 3 and the associated angular momentum redistribution 
at two different times, r ~ 2 (left) and r ~ 6 (right). Note that at earlier times shock 2 
has only one branch (left), whereas it isplits into two at later times (right). Sample gas 
streamlines are shown as black curves with arrows; shocks are shown in red. The large black 
disk is the black hole; the thick black curves without arrows are rough indications of the 
boundaries of the flow. Relative angular momentum is color-coded as shown. The contact 
discontinuity between shocked gas-1 and gas-2 is shown by a dashed black curve. Separations 
are not drawn to scale in order to emphasize the sequence of events. 









When gas-2 passes through this front, it is deflected in a way that causes it to move more 
radially inward; as it does so, it is accelerated by gravity, and its pressure falls by adiabatic 
expansion. When the separation between the shock front and the contact discontinuity grows 
large enough, the shocked gas loses enough pressure by the time it approaches the contact 
discontinuity that it shocks for a second time (right panel of Fig. 11). From this point 
onward, shock 2 is divided into two branches, which we call 2b (the one farther upstream) 
and 2a (the one close to shock 3). 


The first passage through shock 2 of the mass associated with the peak in the mass- 
return rate leads to a surprising effect: the density concentration becomes a quasi-coherent 
object and traces an eccentric orbit around the black hole (its radial oscillation is visible in 
Fig. [I9|. When this “lump” reaches apocenter, the position of shock 3 moves outward; when 
it is near pericenter, shock 3 moves closer to the black hole. 


Angular momentum redistribution is accomplished by deflections at the shocks (Rosswog 


et al. (2009) also report a shock that appears to resemble shock 2 in its early stages and 


remark briefly that it redistributes angular momentum). The shock velocity of shock 2, 
both before and after its split, is directed roughly azimuthally and against the sense of the 
gas orbits; it therefore reduces the angular momentum of gas-2 and deflects it inward. The 
angular momentum lost by gas-2 is transferred to gas-1, deflecting it outward. As a corollary, 
gas-2 also does work on gas-1, reducing its own orbital energy while increasing that of gas-1. 
However, this picture becomes a bit more complicated when shock 2 splits into two separate 
shocks. During this period, angular momentum is transferred from incoming gas-2 at smaller 
radii to gas-2 shocked by shock 2b at larger radii that has then turned inward. Only after a 
second deflection at shock 2a does the angular momentum and energy reach gas-1. Of course, 
gas-1 becomes gas-2 after taking a turn through the pericenter region, so the net result of 
this process is to broaden the distribution of specific angular momentum. Figure 12 shows 
how the box-integrated angular momentum distribution changes over time as a result of 
these processes. From its initial very narrow distribution, the spread in angular momentum 
grows substantially by r = 3. At later times, the growth in width slows, but, as more mass 
at the low end of the angular momentum distribution falls through the inner boundary, the 
average value of the specific angular momentum (vertical lines in the figure) shifts toward 
higher values. The energy distribution also widens, particularly toward more negative orbital 
energy, but because it begins far broader than the angular momentum distribution, the effect 
is less noteworthy. 
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Fig. 12.— Mass distribution of specific angular momentum in the global simulation. Dif¬ 
ferent lines represent the distribution at different times: r = 0.6 (black), 1 (red), 3 (green), 
6 (blue), and 12 (cyan). Broken lines show the average value over the distribution of corre¬ 
sponding color, dm/dj is dimensionless because m is in units of Mbh while j is in units of 

( r /Rg)M c )- 
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4.4. Heating 


To determine the heating rates of the different shocks, we compute the divergence of the 
heat flux carried by the fluid and integrate over a volume containing the shock in question. 
This method suffers from certain inaccuracies (mixing heating by adiabatic compression with 
true entropy generation, by-eye choice of the specific volumes for integration), but in these 
circumstances appears to provide the best feasible option. 


When we apply this method to our data, we find that different shocks create the most 
heat at different times (Fig. 13). Shock 1 forms first, so it dominates the heating early on. 
Once shocks 2 and 3 emerge, their heating rate grows. Shock 3 is a strong shock, with 
Mach number ~ 10, and heats the gas from very low temperatures to ~ 0.1 x the local 
virial temperature. Shock 2 has a lower Mach number because its upstream gas has already 
been heated in shock 1, but nonetheless produces comparable heating. Unfortunately, until 
r ~ 5, shocks 2 and 3 are so close together we cannot reliably separate the volumes affected 
by them in order to compute their heating rates independently. Nonetheless, it is possible 
to say that by r ~ 3-4, the heating in shocks 2 and 3 combined grows to about half that in 
shock 1, and is almost as large as the heating in shock 1 by r ~ 5. 


During this entire earlier period (r < 5), the total heating rate grows, increasing by a 
factor ~ 4 between r ~ 2 and r ~ 5. It remains at about the same level (~ 4 x 10 47 erg/s for 
the parameters of our simulation) until r > 8. However, the dominant contributors change 
after r ~ 5. Heating by shock 1 falls sharply because this is the time at which its Mach 
number becomes ~ 1. Meanwhile, shocks 2 extend to smaller radius, where the larger orbital 
speeds permit their heating rates to grow. This extension also leads to such close proximity 
to the dying shock 1 that we cannot cleanly separate the heating rate in shock 1 from that 
in shocks 2 from r ~ 6 until r ~ 8. During the period after r ~ 3, the rate at which newly- 
arriving matter reaches shock 3 diminishes, while the mass flux passing through the other 
shocks does not decrease nearly as much because matter passes through them repeatedly. 
Consequently, the heating in shock 3, although significant earlier, becomes negligible after 
r ~ 5. Thus, from r ~ 8 onward, the heating is primarily associated with shocks 2. 


As a result of all this heating, the total internal energy grows rapidly until r ~ 5; 
afterward, it is roughly constant (Figure 14). This is because most of the shock heating 
takes place at small enough radii that new heat can be drained through the inner boundary 
by the radial inflow. 


Although the total internal energy generated by the end of the simulation is ~ 300 times 
more than the initial internal energy, it still is only ~ 0.1% of the rest-mass energy of the star 
because the shocks responsible for the heating occur predominantly at large radius. For gas 
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Fig. 13.— Heating rate as a function of time. Total (black curve); shock 1 (red solid curve); 
shock(s) 2 (blue solid curve); shock 3 (solid green curve); shocks 2 and 3 together (dashed 
cyan curve); shocks 1 and 2 together (dashed blue curve). 
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to settle near R p , it would need to lose an order of magnitude more heat. In this sense, the 
hydrodynamic processes taking place in this simulation do not solve the dissipation problem 
posed in the Introduction. It also follows immediately that the heating we have studied 
can account for at most only a very small fraction of the total energy released in a tidal 
disruption event if ultimate accretion onto the black hole has the usual relativistic radiative 
efficiency ~ 10% in rest-mass units. 


4.5. Circularization 


Classical accretion disks are, in most instances, imagined as very thin structures in 
which the material travels on very nearly circular orbits while moving inward on a timescale 
much longer than an orbital period. As a result, departures from axisymmetry are quite 
small. In a tidal disruption event, it is an interesting question whether, and to what degree, 
such a state is achieved. To measure progress toward such “circularization”, we define several 
independent measures. 


The first is the level of nonaxisymmetry in its surface density. Figure [15] shows the 
rms fractional surface density fluctuation around the azimuthal coordinate at each radius, 
<5(r) = [£(r, </>)/S(r) — l] rms , at several times. Here S is the mean over azimuthal angle qb. At 
large radii (r > 400M-600M), 5 remains ~ 1 until very nearly the end of the simulation; this 
large asymmetry reflects the narrowness of the tidal streams. The large asymmetry at all 
radii when r = 1 is due to the same cause. Closer in (100M < r < 400M), the disk remains 
quite asymmetric (5 ~ 0.5) until r > 4, when 6 falls to ~ 0.2-0.3 and stays at that level 
until the end of the simulation. In this region, the asymmetries are largely associated with 
shocks 1, 2a, and 2b. At still smaller radii (30M < r < 50M), there is a modest increase 
in asymmetry that arises from the inward tunneling of particularly low angular momentum 
gas. The upshot from this measure is that the flow is quite far from circular at least until 
r ~ 4, and remains significantly non-circular even at r ~ 12. 


Our second measure is the ratio of radial to azimuthal velocity Ity/v^. Whereas 5(r) 
measures morphological circularity, \v r /v^,\ measures the circularity of orbital motions. Fig¬ 


ure 


16 shows the radial profile of the density-weighted average of Ity/i^l, 


Rv{r) = [Vr/v^, 


J p^r/v^rdepdz 
f prd(j)dz 


( 7 ) 


As can be seen in this figure, the degree of orbital circularity improves steadily over time. 
Early on (r ~ 1), R v has a rather complex dependence on radius; this dependence is due 
to the initial tidal stream structure. From r ~ 2 onward, R v declines steadily at all radii 
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Fig. 14.— Time evolution of total internal energy inside simulation box (green), total 
internal energy escaped from the simulation box (red), and sum of them (black). 
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Fig. 15.— Radial profile of rms azimuthal fractional surface density fluctuation, S(r) = 
[S(r, 0)/E(r) — l] rms at several times. 
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from ~ 1 to ~ 0.2-0.4 by the end of the simulation. Just as steadily, the radius at which 
R v is minimum moves outward, reflecting the progressively larger portion of the flow that is 
at least partially circularized. In fact, judged on an absolute standard, the fluid orbits are 
approaching a reasonable degree of circularity by late times. The relative magnitude of the 
radial velocity would also be roughly in line with expectations for the mean inflow speed if a 
disk this thick were subject to the MHD stresses that ordinarily drive accretion, even though 
these stresses are absent here. 


A third measure is given by the nominal eccentricity of the fluids’ orbits; that is, the 
eccentricity of a test-particle orbit with the same specific angular momentum and orbital 
energy: 

e=[l + 2 ej 2 l(GM) 2 )] 1/2 , (8) 

where the specific orbital energy e = (l/2)u 2 — GM/r and the specific angular momentum 
j = rVf/). For our diagnostic, we construct its mass-weighted azimuthal and vertical average 


e r (r) = 


f perdcfxiz 
f prdcfidz 


(9) 


Figure 17 shows the radial profile of e r for several times ranging from the very beginning 
of the simulation to the end. This measure also decreases steadily over the duration of the 
simulation, although its rate of decrease decelerates after r ~ 6. By the end of the simulation, 
the mean eccentricity has diminished to ~ 0.3-0.4 over most of the flow’s volume. In other 
words, this diagnostic shows a level of circularity similar to the one implied by R v . 


Lastly, we examine the proportion of rotational support against gravity as measured by 
the ratio between the gas’s specific angular momentum and that required for a circular orbit 
with that fluid element’s instantaneous radius; in a classical disk, this ratio would be unity. 
Figure 18 shows the mass-weighted average of this quantity as a function of time at several 
radii. 

f pjrd(j)dz 


J(r) = jmean(r)/j K (r) = 


ji<{r) f prd(j)dz' 


( 10 ) 


Initially, almost complete rotational support is implied at radii near and inside r p (60 M < 
r < 100M), while rotational support plays a much smaller role everywhere else. However, 
over time all locations converge toward a value of J ~ 0.8, a bit greater for r < 100M, 
somewhat less at r > 1000M. This much angular momentum is enough to create a somewhat 
flattened structure, but 1 —J is large enough that substantial pressure support is also implied. 
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Fig. 16.— Radial profile of the mass-weighted azimuthal average of the ratio of radial to 
azimuthal velocity. Black, red, green, blue, and cyan colors correspond to times r =1, 3, 6, 
9, and 12, respectively. 
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Fig. 17.— Radial profile of mass-weighted mean eccentricity. Black, red, green, blue, and 
cyan colors correspond to r =1, 3, 6, 9, and 12, respectively. 
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Fig. 18.— Radial profile of density weighted average of J, the ratio of the mass-weighted 
specific angular momentum to that required for a circular orbit, at several times. Black, red, 
green, blue, and cyan colors correspond to r =1, 3, 6, 9, and 12, respectively. 
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4.6. Actual mass accumulation rate 


As discussed in Sec. [lj it has long been assumed that in tidal disruption events the 
ballistic mass-return rate (after a quick rise to a maximum at r ~ 1) is well approximated 
by (M*/3 )t - 5 / 3 , and that this same expression is a good predictor of both the rate Mdisk at 
which matter accumulates in an accretion disk of radius ~ 2 R p and the rate M acc at which 
matter enters the central black hole. We have already shown (Sec. 3.1) that the ballistic 
mass-return rate reaches its maximum at r ~ 1.5 and asymptotes to a curve oc r~ 5 / 3 , but 
its maximum value is about twice the value predicted by the r~ 5 / 3 curve at the time of that 
maximum and asymptotes to r _5//3 only for r > 3. 


Our hydrodynamic results indicate stronger contrasts with the usual assumptions about 
Mdisk and M acc . First of all, the “disk” structure forms at considerably larger radius than 


~ 2 Rp. As shown in Figure 19, even at r ~ 12, only ~ 1/3 of the bound mass has moved 


within 2 R p of the black hole. In fact, to reach the point within which half the bound mass 
can be found requires moving out to r ~ 6 R p , and only ~ 3/4 of the bound mass is within 
16R p at r = 12. Thus, most of the mass accumulates at scales comparable to the semi-major 
axis of the most bound material a m ; n (in this case, ~ 500M = 5 R p ), rather than near r ~ R p . 

The extended mass distribution and its azimuthal asymmetry together make the accu¬ 
mulation of mass in the accretion flow considerably slower than the mass-return rate would 
indicate. Half the ultimate mass within 2 R p arrives after r ~ 5; even bringing half the 
ultimate mass within 8 R p isn’t accomplished until r ~ 3. Second, and more surprisingly, 
M(< r) oscillates as a function of time, particularly for 2 R p < r < 10 R pi i.e., Mdisk changes 


sign frequently across that entire range of radii (Fig. 20). At early times, a great deal of 


mass moves within 10 R p as the tidal streams begin their approach; however, after passing 
near pericenter, the matter in those streams swing back out, reaching as far as 10_R P (see 
also the first panel of Fig. [6|. This effect leads to an oscillation in the mass within ~ 8 R p 
because the streams are not evenly spread in orbital phase—there is a concentration in the 
lump mentioned earlier. Thus, any measure of Mdisk follows the ballistic return-rate curve 
only with sizable fluctuations, both above and below the ballistic return-rate, until r > 5. 


Stream deflection of shocks does, however, partially aid in directing mass inward, sup¬ 
plementing the usual source of internal stress in accretion flows, correlated MHD turbulence. 
However, the rate at which this mechanism drives matter through our inward boundary is 
well below the ballistic return rate until r > 5, and the rate reached at that time is only 
~ O.lx the peak rate of mass-return as predicted by a purely ballistic model. This rate, 
of course, is still not the rate M acc at which mass will flow into the black hole; that rate 
depends on the strength of further angular momentum transport, presumably due to MHD 
turbulence. 





M <r [MJ 
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Fig. 19.— Mass interior to a given radius as a function of time. The various radii plotted 
are: 16_R P (yellow curve), 8 R p (magenta curve), 4 R p (cyan curve), 2 R p (red curve), R p (green 
curve), r in (blue curve). 
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These results contrast strongly with classical expectations. By the end of the simulation 
at r ~ 12, 28% of the bound mass has passed through the inner boundary. In the conventional 
view, in which the accretion rate matches the return rate’s %/ 3 time-dependence, this much 
mass would have been accreted by r ~ 1.6, and 80% would have been accreted by r ~ 12. 


4.7. Cooling time 

One of the key assumptions of our calculation was that the flow is adiabatic. We can 
test the validity of that assumption by estimating the radiative cooling time f coo i rs_/ t t H/c, 
where tt is the (half-)optical depth to Thomson scattering in the vertical direction. The 
cooling time declines outward, but rather slowly, because H/r varies little with radius. As 
a result, the outward decline in surface density is partially offset by the outward increase in 
r/c. For our parameters, t coo i is ~ 50-150 yr, or ~ 5 x 10%, for r < 1000M. 


5. Discussion 

5.1. How does the accretion rate onto the black hole vary in time? 

As we have shown, shock deflection removes angular momentum and orbital energy from 
parts of the stellar debris, transferring this angular momentum and energy to other parts. 
By this process, roughly a third of the bound mass is forced inward relatively quickly, albeit 
on a timescale rather slower than the nominal mass-return timescale. Nearly all the mass 
that finds its way to r < 2% has, by the end of the simulation, gone all the way through 
the inner boundary at 30% = 0.3%. However, half the inflow through the inner boundary 
occurs after r ~ 6. Once within 30%, if MHD turbulence develops in the usual way, the 
inflow time is relatively swift: 

inflow ~7x 10 3 (a/0.1)- 1 [(%r)/0.5]- 2 (r/30%) 3 / 2 %/c (11) 

or Tinflow ~ 0.1(a/0.1) -1 [(%/r)/0.5] -2 (r/30%) 3 / 2 . Here a is the usual ratio of vertically- 
integrated stress to vertically-integrated pressure. Because tinflow of the matter pushed early 
through the inner boundary is much shorter than the time it takes to reach that boundary, 
accretion of this material onto the black hole and arrival at the inner boundary are nearly 
instantaneous. 

The majority of the bound mass, meanwhile, settles into the flow at much larger distance. 
At the characteristic radius of this portion of the stellar debris, the inflow time is much longer: 

^inflow ~ 1-3 x 10 6 (a/0.1)- 1 [(%r)/0.5]- 2 (r/1000%) 3 / 2 %/c, (12) 
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Fig. 20.— Several different measures of the mass inflow rate: ballistic, i.e., pseudo-particles 
(black curve); the mass flow through the inner boundary (M acc , blue curve); the rate of 
increase of the mass within 2 R p (green curve); and the rate of increase of the mass within 
10R p (cyan curve). There are frequent episodes of net outflow, these are denoted by dashed 
curves whose color corresponds to the radius of reference. 
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or Tinflow ~ 19(a/0.1) 1 [(H/r)/ 0.5] 2 (r/1000i? 3 ) 3 / 2 . Thus, the majority of the mass arrives 
at the black hole only after a considerable delay. 


We can roughly predict the accretion rate at later times by constructing a crude ap¬ 
proximate description for the action of turbulent MHD stresses. Suppose that the time at 
which matter accumulated at radius r reaches the black hole’s IS CO is 


4cc = 4rrive + [« 2 (H / r) 2 + 7l sat ] (r / R g f /2 R g / C. 
The rate at which matter accretes onto the black hole is then roughly 

dM ,, dr 

= 2nrY,[r(t acc )]- 


dt dt^ 

where we have identified the time of observation t with 4 


(13) 


(14) 


Several parameters govern the behavior of this model. Our hydrodynamic simulation 
results suggest that the time at which matter settles into the flow, f arr i ve , is — 12 1 0 . Although, 
as shown by Figure 19, M(< r ) approaches its final value at most radii by r ~ 5, the data 
displayed in Figure [6] demonstrate that the motion of the matter remains dominated by 
shocks and non-circular flow until nearly the end of the simulation. Because this is a newly- 
formed flow, a finite time is required for the MHD turbulence to reach saturation; we call 
this n sat dynamical times. Numerous simulations, both of shearing boxes and global disks, 
have shown that ~ 10 orbital periods are generally required for the turbulence to reach 
saturation; on this basis, we estimate n sa t ~ 50. In a steady-state disk, the mean inflow time 
is a~ 2 (H/r)~ 2 (r / R g ) 3 / 2 R g /c. In this context, the saturation level of the MHD turbulence, 
and therefore the best value of a for estimates, is more uncertain than usual. Even at 
r > 10, the accretion flow is still not entirely circular or axisymmetric. More importantly, 
it is intrinsically time-dependent. We therefore suppose only that a ~ 0.01-0.1. 


The predicted accretion rate derived from such a model is shown in Figure |2l| The peak 
accretion rate is due to the early forcing of matter inward. It is reached rather later (r ~ 3 
rather than r ~ 1) than has generally been assumed and is maintained rather longer (until 
r ~ 8). The value of the peak accretion rate is ~ O.lx the peak in the ballistic mass-return 
rate, the usual estimator of the accretion rate. This order of magnitude reduction in the 
maximum mass-return rate is due in part to the fact that the width of the peak is greater 
by a factor ~ 5 and in part to the fact that the mass arriving during the peak is only about 
1/2 the mass delivered at the peak ballistically. After r ~ 8, accretion of the early-arriving 
material drops off very quickly; by r ~ 10-15 (the exact value depending on a(H/r) 2 ), 
accretion of the material initially placed at larger distance dominates. The accretion rate 
at late times also depends on a(H/r) 2 , but is likely to be at least a factor of several below 
the peak accretion rate resulting from the early accretion. On the other hand, once the 
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later mechanism becomes established, the accretion rate can remain roughly constant for a 
significant time, ~ 30-100 in r units (multiples of t 0 , the orbital period of the most tightly- 
bound tidal stream). At truly long times, the accretion rate resumes a power-law drop-off 
whose slope depends on the radial surface density profile at the time the mass settles into 
the flow and on the stress parameters. 


If the accretion flow is radiatively efficient, the lightcurve should closely mirror the ac¬ 
cretion rate. It is convenient in this context to translate the dimensionless units of Figure [2l] 
into the dimensionless units of ratio to the Eddington accretion rate Me- the conversion 
factor is such that 10 _7 M*/(GM BH /c 3 ) = 5.57 x 10 7 (?7/0.1)M£, where r] is the radiative effi¬ 
ciency of accretion in rest-mass units. Thus, for the parameters of our simulation, the peak 
accretion rate is extremely super-Eddington, ~ 1 x 10 8 Me- For such a very large accretion 


rate, photon trapping effects may severely limit the flow’s radiative efficiency (Begelman 


1979 Abramowicz et al. 1988). On the other hand, magnetic buoyancy effects (Blaes et al. 


2011) may be able to carry heat to the outer layers rapidly enough to restore some or perhaps 


all of the system’s radiative efficiency (Jiang et al. 2014). 


Despite these uncertainties, it is worthwhile considering what our approximate model for 
the accretion rate predicts in the event of efficient radiation. Particularly if accretion stresses 


are relatively large (e.g., the green curve in Fig. 21), the form of the lightcurve remains one in 


which there is a rapid rise to a peak followed by a long, roughly power-law decay. The shape 
of this decay is, however, not exactly the conventional f -5 / 3 decline. Although there are 
periods when it is reasonably well described by a power-law with approximately this index, 
there is also a period during which the accretion rate is roughly constant at a level that is a 
fraction of the peak value. The internal accretion rate, parameterized by a(H/r) 2 , plays the 
greatest role in determining the level at which this flat section occurs and its duration. To 
the extent that photon-trapping plays a role, the peak of the lightcurve would be depressed 
and the early phase of the lightcurve flattened. 


5.2. Extrapolation to main sequence star disruptions and larger A/bh 

The parameter regime of our simulation—a white dwarf passing close to a 5OOM 0 black 
hole—refers to a minority subclass of possible tidal disruption events. It is worthwhile 
examining the extent to which its implications can be extrapolated to other, more common, 
parameters. The most natural way to make the connection is through R p in gravitational 
units. We studied an event in which a white dwarf has R p = Rt = 100i? 3 ; if a main sequence 
star of mass Mms were to have an encounter in which this were also true, the black hole 
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Fig. 21.— Accretion rate onto the black hole. The black curve shows the usual assumption, 
in which the accretion rate is the same as the ballistic return rate. The red and green 
curves show the sum of the rate at which matter passes through the inner boundary in our 
simulation and the accretion rate predicted by our approximate model. They differ in their 
choices of a, 0.1 for the green curve, 0.01 for the red curve. 
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mass would be 

M' b h = M bh (M wd /M ms ) 1/2 (R ms /R wd ) 3/2 . (15) 

Because main sequence star radii are ~ 100 x larger than white dwarf radii, the black hole 
for the corresponding main sequence tidal disruption would be ~ 1000 x larger than the 
black hole in our simulation; to be specific, Mg H = 3 x 10 5 M M g because i£* oc M* for stars 
with mass ~ M 0 on the main sequence. 

Although we have matched R p /R g , the orbital scale for the most tightly-bound tidal 
streams will not be the same: it is a m i n /R g ~ (M^/M^iRt/Rg). Thus, in gravitational 
units, a min for a main sequence disruption is an order of magnitude larger than for a white 
dwarf disruption. At early times, the speeds of shocks 2 ~ (a min /Rg)^ 1 ^ 2 , so the heat 
generated per unit mass passing though the shocks would fall by about that same single order 
of magnitude, although as they extend inward at later times, this effect will be somewhat 
mitigated. 

On the other hand, for fixed R p /R g , the photon diffusion time in shocked material scales 
tx a m } n in physical units, or oc M BH ~ 4//3 . We would then expect a decrease in the cooling time 
by a factor ~ 10 4 , to ~ 3 x 10 5 s. The relevant standard of comparison for the cooling time 
is t 0 , and it scales oc M^ H 3 ^ 2 for fixed R p /R g . Consequently, the ratio t coo i/t 0 is extremely 
sensitive to black hole mass, oc Mg H ~ ll,/6 . For this reason, the adiabatic approximation, 
although well-supported for our parameters, breaks down for Mg H > 2 x 1O 4 M 0 . Post¬ 
shock cooling when the black hole mass is larger should then be significant, possibly further 
postponing the erasure of density inhomogeneities in the flow. Quick post-shock cooling 
could also create a precursor signal preceding the main flare that takes place when matter 
begins to accrete onto the black hole. However, the available energy per unit mass is at 
most ~ c 2 R g /a m in , and this will always be small compared to the black hole accretion energy 
release unless nearly all of the photons created near the black hole are trapped in the inflow. 


Our most important result is impervious to these changes of scale. The orbital period of 
the most tightly-bound matter t 0 determines all aspects of the inflow: the ballistic dynamics 
of the tidal streams, the orbital period at the characteristic radius of shocks 2 and 3, and 
the inflow time from that characteristic radius (tmfiow ~ oi~ l (H / r)~ 2 to). Because the mass 
of the disrupted star is ~ M 0 whether it is a white dwarf or on the main sequence, the 
time-dependence of the accretion rate depends only on this single parameter. Therefore, by 
plotting the accretion rate in terms of time in units of to, the time-dependence shown in 
Figure 21 should apply reasonably well to all TDEs, with only minor adjustments due to 
variations in R p /R t , the internal structure of the disrupted star, etc. 


Thus, all that is necessary to predict the time-dependence of accretion rate during a 
main sequence TDE is to reinterpret r in terms of the appropriate t 0 . Comparing our results 
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to the predictions of the classical theory, we find that we, too, expect a sharp rise followed 
by an extended, roughly power-law decay. However, there are two contrasts. For a given 
set of parameters (R p /R tj Mbh? M*), the peak accretion rate is ~ lOx smaller, and the 
relation between the duration of the peak is about 5x longer. The factor of ten diminution 
in maximum accretion rate means that it becomes ~ 80M* ^bhg i n Eddington units. 
This reduction in the accretion rate relative to Eddington should also lead to a weakening of 
photon-trapping. Because t 0 oc (M^Ml ) b 2 , the geometric mean of the black hole and 
stellar masses inferred from a measured timescale is ~ 5 x smaller than the one predicted by 
the classical theory. 


6. Summary 


When a star is tidally disrupted by a black hole, its matter is dispersed into thin streams 
that travel far from the place where the star was torn apart before returning to the vicinity 
of the black hole. Before those streams can accrete onto the black hole and generate the 
energy for a flare, their orbital energy must be dissipated, much of the heat generated must 
be lost, and their angular momentum must be reduced to less than that supporting an orbit 
at the black hole’s ISCO. In this paper, we have focused on the first step—how the orbital 
energy is dissipated—and we have done so by consciously ignoring MHD effects, which almost 
certainly are important to the last step (angular momentum transport) and possibly even 
to the middle step (heat loss: see Blaes et al. (201 1[) ; Jiang et ah] (2014)). 


Previous work has emphasized the role of what we have named “shock 1” (also called 
the “nozzle” shock), a shock formed near the pericenter of the star’s orbit, where the various 
tidal streams converge. We have identified a shock system, shocks 2 and 3, that appears near 
the apocenter of the streams’ orbits. These shocks are created by an intersection between 
a stream that has just swung through pericenter and another stream that is falling toward 
its pericenter. If the orbits of the two streams have a modest (~ 10°) contrast in apsidal 
angle, the angle between their velocities at the point where their orbits intersect can be much 
larger, ~ 90°. We have shown that just such a swing in apsidal angle is created by relativistic 
effects during the disruption, particularly the relativistic precession in apsidal angle of the 
star’s trajectory. Because R t ~ O(10)i? 9 is typical, these modest contrasts in apsidal angle 
can be expected generically in stellar tidal disruptions. 


All three shocks dissipate orbital energy, although the amount they dissipate is small in 
rest-mass terms, ~ 5 x 10 50 erg altogether in this simulation, ~ 10 _3 x the rest-mass energy 
of the stellar mass bound to the black hole. Although shocks 2ab contribute the largest 
share of the heating, shock l’s share is smaller by only a factor of a few, and part of this 







-48 - 

shock pair moves in toward the location of shock 1 during the time that their heating rate 
is greatest. We therefore expect the amount of heating in these shocks to scale with R g /R p , 
a number that is relatively insensitive to black hole mass. 

In addition, by deflecting the tidal streams, the shocks substantially widen what was 
initially a rather narrow specific angular momentum distribution, causing some material to 
plunge more directly toward the black hole while other material follows a more circular orbit 
at larger distance. This mechanism can cause an interesting fraction of the mass to reach 
radii a factor of a few smaller than R p without any of the usual MHD internal stresses. The 
broader angular momentum distribution also leads to a wider range of orbital paths for gas 
passing through the pericenter region for a second or third time, significantly broadening the 
front associated with shock 1, extending it both inward and outward radially. 

Previous work had also assumed that the returning tidal streams would quickly join 
an accretion disk whose outer edge would lie at the radius, 2 R p , where the specific angular 
momentum to support a circular orbit matched the specific angular momentum of the star’s 
trajectory. However, we find that shocks 2 and 3 lead to such a build-up of material in 
the vicinity of r ~ a m j n , the semi-major axis of the most tightly-bound material, that the 
outer edge of the accumulated material is on that scale, rather than the considerably smaller 
Rp. Moreover, there are significant structural contrasts between this “disk” and a classical 
accretion disk. Even long after the great majority of the star’s bound mass has returned to 
pericenter at least once, the structure is strongly asymmetric, and the fluid orbits are rather 
eccentric. 

These results have strong implications for the interpretation of TDE observations. Ar¬ 
rival of matter at radii close enough to the black hole to release a substantial amount of 
energy appears to be subject to several kinds of delays. A significant minority of the star’s 
bound mass is pushed well within the star’s pericenter, but only after a time ~ 5fo, where to 
is the orbital period of the most tightly-bound matter, whereas classical expectations about 
tidal disruptions assumed that much of the star’s mass would reach the black hole only 1- 
21 0 after the star passed through pericenter. Accretion of the majority of the star’s bound 
mass is delayed by a considerably longer time because, at its scale ~ a min , inflow driven by 
the usual internal stresses associated with MHD turbulence takes a much longer time. Al¬ 
though the classical scaling of characteristic timescale oc (AdBuM*) 1 ^ 2 remains appropriate, 
these several sources of delay insert a proportionality constant in this relation ~ 5x larger 
than generally estimated. As a result, the product of Mbh and M* inferred from a particular 
lightcurve would be an overestimate by ~ 25, and the peak accretion rate in Eddington units 
is reduced by a factor ~ 10. However, this suggestion requires confirmation (or refutation) 
by a future MHD simulation. It is conceivable that MHD effects (e.g., stresses associated 
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with magnetic field amplification by fluid shear) could also influence the flow even before 
formation of anything resembling an accretion disk. 
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